


FUNCTIONAL DETERMINANTS IN HIGHER DIMENSIONS 
USING CONTOUR INTEGRALS 

KLAUS KIRSTEN 

o 

^vj , Abstract. In this contribution wc first summarize how contour integration 

methods can be used to derive closed formulae for functional determinants 
of ordinary differential operators. Wo then generalize our considerations to 
partial differential operators. Examples are used to show that also in higher 
dimensions closed answers can be obtained as long as the eigenvalues of the 
differential operators are determined by transcendental equations. Examples 
considered comprise of the finite temperature Casimir effect on a ball and the 
functional determinant of the Laplacian on a two-dimensional torus. 



p i' 1. Introduction 

D : 

(~| ^ Functional determinants of second order differential operators are of great im- 

portance in many different fields. In physics, functional determinants provide the 
so-called one-loop approximation to quantum field theories in the path integral for- 

^ . mulation [211 El] • In mathematics it describes the analytical torsion of a manifold 

in: EH- 

0^ ■ Although there are various ways to evaluate functional determinants, the zeta 

function scheme seems to be the most elegant technique to use [Qj [HI [TTJ [31] . 
This is the method introduced by Ray and Singer to define analytical torsion [17] . 

^O I In physics its origin goes back to ambiguities in dimensional regularization when 

applied to quantum field theory in curved spacetime [TTJ [29] . 

For many second order ordinary differential operators surprisingly simple answers 
can be given. The determinants for these situations have been related to boundary 
values of solutions of the operators, see, e.g., PfT0 l [T2 ] [22 l [23 l [26 l [36 l [39 l l40]. 
Recently, these results have been rederived with a simple and accessible method 
?— ( ' which uses contour integration techniques [33l [34l [35] . The main advantage of this 

approach is that it can be easily applied to general kinds of boundary conditions 
[35] and also to cases where the operator has zero modes [34l [35]; see also [371 
l38l [42]. Equally important, for some higher dimensional situations the task of 
finding functional determinants remains feasible. Once again closed answers can be 
found but compared to one dimension technicalities are significantly more involved 
[131 [T4] . It is the aim of this article to choose specific higher dimensional examples 
where technical problems remain somewhat confined. The intention is to illustrate 
that also for higher dimensional situations closed answers can be obtained which 
are easily evaluated numerically. 

The outline of this paper is as follows. In Section 2 the essential ideas are pre- 
sented for ordinary differential operators. In Section 3 and 4 examples of functional 
determinants for partial differential operators are considered. The determinant in 
Section 3 describes the finite temperature Casimir effect of a massive scalar field 
in the presence of a spherical shell [24l [25] . The calculation in Section 4 describes 
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determinants for strings on world-sheets that are tori [46l [50] and it gives an aher- 
native derivation of known resuhs. Section 5 summarizes the main results. 

2. Contour integral formulation of zeta functions 

In this section we review the basic ideas that lead to a suitable contour inte- 
gral representation of zeta functions associated with ordinary differential operators. 
This will form the basis of the considerations for partial differential operators to 
follow later. 

Wc consider the simple class of differential operators 

on the interval / = [0,1], where V{x) is a smooth potential. For simplicity we 
consider Dirichlet boundary conditions. From spectral theory |41j it is known that 
there is a spectral resolution {(/)„, Xn}^=i satisfying 

PM^) = KMx), MO) = 0n(i) = 0. 

The spectral zeta function associated with this problem is then defined by 
"^^^ (2.1) Cp(s) = ^A-^ 

n=l 

where by Weyl's theorem about the asymptotic behavior of eigenvalues [35] this 
series is known to converge for 5fts > 1/2. 

If the potential is not a very simple one, eigenfunctions and eigenvalues will 
not be known explicitly. So how can the zeta function in equation (j2.ip . and in 
particular the determinant of P defined via 

det P = e-'^'^^°\ 

be analyzed? From complex analysis it is known that series can often be evaluated 
with the help of the argument principle or Cauchy's residue theorem by rewriting 
them as contour integrals. In the given context this can be achieved as follows. Let 
A e C be an arbitrary complex number. From the theory of ordinary differential 
equations it is known that the initial value problem 

resOlal (2.2) (F - X)ux{x) = 0, Ux{0) = 0, <(0) = 1, 

has a unique solution. The connection with the boundary value problem is made 
by observing that the eigenvalues A„ follow as solutions to the equation 

res02 | (2.3) ua(1) = 0; 

note that u\{l) is an analytic function of A. 

With the help of the argument principle, equation (|2.3p can be used to write the 
zeta function, equation (|2.ip . as 



"^^^ (2.4) ^p(s)= /rfAA-^— InuA(l). 

' ' z-Ki J d\ 



Here, 7 is a counterclockwise contour and encloses all eigenvalues which we assume 
to be positive, see Figure 1. The pertinent remarks when finitely many eigenvalues 
are non-positive are given in |35j . 
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A-plane 



Figure 1. Contour 7 used in equation (|2.4p . 



The asymptotic behavior of liA(l) as |A| — > 00. namely 



«a(1) 



1\/A 



VA 



imphes that this representation is vahd for 5Rs > 1/2. In order to find the determi- 
nant of P we need to construct the analytical continuation of equation ()2.4p to a 
neighborhood about s = 0. This is best done by deforming the contour to enclose 
the branch cut along the negative real axis and then shrinking it to the negative 
real axis, see Figure 2. 



A-plane 



Figure 2. Contour 7 used in equation (|2.4[) after deformation. 
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The outcome is 

oo 

Tiio^ (2.5) Cp{s)^ / rfA A-^— lnu_A(l). 

TT J d\ 



To see where this representation is well defined notice that for A — > oo the behavior 
follows from [IT] 

The integrand, to leading order in A, therefore behaves like A * ^/^ and conver- 
gence at infinity is established for 5fts > 1/2. As A — >■ the behavior A~^ follows. 
Therefore, in summary, (|2.5p is well defined for 1/2 < 5fts < 1. In order to shift the 
range of convergence to the left we add and subtract the leading A — > oo asymptotic 
behavior of u_a(1). The whole point of this procedure will be to obtain one piece 
that at s = is finite, and another piece for which the analytical continuation can 
be easily constructed. 

Given we want to improve the A — > oo behavior without worsening the A — )■ 
behavior, we split the integration range. In detail we write 

res05 | (2.6) Cp(s) = Cpj{s) + Cp^asis), 



where 
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Cp,as{s) 


= 


sin TTS /■ 


/ 



rfAA-^-^ln«_;,(l) 
dX 



(2.7) + !i^|dAA-^in(,_,(i)2VAe-v^ 



res07 (2.8) (Ip^asis) = d\\-'—ln 



d ' -^ 



e' 



1 



dX \2VX 



By construction, Cpj{^) is analytic about s = and its derivative at s = is 
trivially obtained, 

(2.9) Cp,/(0)=lnu-i(l)-lnuo(l)-ln(u-i(l)2e-i) = -ln('^^^ j . 

Although the representation (|2.8p is only defined for 5Rs > 1/2, the analytic con- 
tinuation to a meromorphic function on the complex plane is found using 

dX X-" = -^— for ^a > 1. 
a- 1 

1 



This shows 



and furthermore 



sin TTS / 1 1 

Cp,as[S) = -^ - - 

2n Vs - 5 s 



Cp,asiO) 
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Adding up, the final answer reads 

resOSa | (2.10) Cp(0) = - ln(2uo(l)). 

For the numerical evaluation of the determinant, not even one eigenvalue is needed. 
The only relevant information is the boundary value of the unique solution to the 
initial value problem 

(—^ + V{x))uo{x)=0, uo{0)=0, u'o{0)^l. 

General boundary conditions can be dealt with as easily. The best formulation 
results by rewriting the second order differential equation as a first order system 
in the usual way. Namely, we define vx{x) = du\{x) / dx such that the differential 
equation (j2.2p turns into 

d fux{x)\ _ f 1 \ fux{xy 



res09 (2.11) -— , / x I - l t// n \ a I 1 / \ 
dx\vx{x)J \ V[x)-X J \vx{x) 

Linear boundary conditions are given in the form 

where M and N are 2x2 matrices whose entries characterize the nature of the 
boundary conditions. For example, the previously described Dirichlet boundary 
conditions are obtained by choosing 

In order to find an implicit equation for the eigenvalues like equation (|2.3|) we use 
the fundamental matrix of p. lip . Let u\^ (x) and m^ (x) be linearly independent 
solutions of (|2.11l) . Suitably normalized, these define the fundamental matrix 



(1)/ \ (2),' \ 

l)/- \ (2)/ \ 



Hx{x)={ %" "f^,:. h iJA(0)=Id2x2. 



The solution of (|2.1ip with initial value {ux{0),vx{0)) is then obtained as 

vx[x)J \vx(0) 

The boundary conditions (|2.12p can therefore be rewritten as 

uxiO)\ /O^ 



|res01l| (2.13) (M + NHxil)),^^^^^. ,p 

This shows that the condition for eigenvalues to exist is 

det{M + NHx{l)) =0, 

which replaces (j2.3p in case of general boundary conditions. The zeta function 
associated with the boundary condition (|2.12p therefore takes the form 

(pis) - — / dA A-^— lndet(A/ + NHx{l)) 
2m J d\ 

7 

and the analysis proceeds from here depending on M and N . HP represents a 
system of operators one can proceed along the same lines. Note, that we have 
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replaced the task of evaluating the determinant of a differential operator by one of 
computing the determinant of a finite matrix. 

The procedure just outlined is by no means confined to be applied to ordinary 
differential operators only. In fact, the zcta function associated with many boundary 
value problems allowing for a separation of variables can be analyzed using this 
contour integral technique. In more detail, starting off with some coordinate system, 
see, e.g., [43], eigenvalues are often determined by 

where j is a suitable quantum number depending on the coordinate system consid- 
ered and Fj is a given special function depending on the coordinate system; e.g. 
for ellipsoidal coordinate systems the relevant special function is the Mathieu func- 
tion. Continuing along the lines described above, denoting by dj an appropriate 
degeneracy that might be present, we write somewhat symbolically 



[^ (2.14) CpW = ^d,^ IdAA-^^ In i^, (A), 



the task being to construct the analytical continuation of this object to s = 0. 
The details of the procedure will depend very much on the properties of the special 
function Fj that enters. For example, on balls Bessel functions are relevant [H|6l[7], 
the spherical suspension [3] , or sphere-disc configurations [27l [32] , involve Legendre 
functions, ellipsoidal boundaries involve Mathieu functions etc. For many examples 
relevant properties of Fj (A) are not available in the literature and need to be derived 
using techniques of asymptotic analysis [lU [44l [45] . For quite common coordinate 
systems like the polar coordinates this is not necessary. When the asymptotics is 
known, the relevant integrals resulting in (|2.14p need to be evaluated and closed 
expressions representing the determinant of partial differential operators are found. 
Although the remaining sums in general cannot be explicitly performed, the results 
obtained are very suitable for numerical evaluation. 

3. Finite temperature Casimir energy on the ball 

Let us now apply the above remarks about higher dimensions using the general 
formalism described in fl4| . As a concrete example we consider the finite temper- 
ature theory of a massive scalar field on the three dimensional ball. Using the zeta 
function scheme we have to consider the eigenvalue problem 

resl| (3.1) P<f>x{T,x):= (--^-A + mAMr,x)^X^MT,x), 

where r is the imaginary time and x E B^ := {x G R-^llxl < 1}. We have written A^ 
for the eigenvalues to avoid the occurrence of square roots in arguments of Bessel 
functions later on. 

For finite temperature theory we impose periodic boundary conditions in the 
imaginary time, 

(J>x{t,x) = 0a(t + /3,x), 

where /3 is the inverse temperature, and for simplicity we choose Dirichlet boundary 
conditions on the boundary of the ball, 

</'A(r,x)||-|^i = 0. 



res2 



lEH (3-3) ^--ol^Cp/^^iO), 
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The zeta function associate with this boundary value problem is then 
(3.2) Cp(s) = E^^'^ 

A 

and the energy of the system is defined by 

ld_ 

2d/3' 

where ^ is an arbitrary parameter with dimension of a mass introduced in order 
to get the correct dimension for the energy. For a full discussion of its relevance 
in the renormalization process in this model at zero temperature see [5]. That 
discussion remains completely unchanged at finite temperature and we will put 
yU = 1 henceforth. 

Given the radial symmetry of the problem we separate variables in polar coor- 
dinates according to 

1 . 2xnx 

with the spherical surface harmonics Y£m{(^, f) [10] solving 

sm 9 oip'' sm of) oO 

and with the Bessel function J^{z), which is the regular solution of the differential 
equation [55] 

Imposing the boundary condition on the unit sphere, 

res4a (3.4) >/f + i (w£ j ) = , 

determines the eigenvalues. Namely, 

/27rn\ ^ 
res4| (3.5) Xl^^^{—-\ +ijjj-+m^, nel., ^ S No, j G N. 

This leads to the analysis of the zeta function 

OO OO OO 

f^^ (3.6) Cp{s)= Y. EE(2^+l)(Pn + ^'.+"^')"'' 

n— — OO ^^0 j — 1 

where we have used the standard abbreviation p„ = {27Tn)/f3. The factor {21 + 1) 
represents the multiplicity of eigenvalues for angular momentum i. 

The zeroes coij of the Bessel functions Jg_^_i{u}ij) are not known in closed form 
and thus we represent the j-summation using contour integrals. Starting with 
equation p.4p and following the argumentation of the previous section, this gives 
the identity 

OO OO r ^ \ ^ 

[^^^ (3.7) Cp(s)= E E(2^+1)/ ^(P» + ^' + "^'r'dA^'''^^+^^^)' 

valid for Sfts > 2. The contour 7 runs counterclockwise and must enclose all the 
solutions of p.4p on the positive real axis. The next step is to shift the contour 
and place it along the imaginary axis. As A — >■ we observe that to leading order 
J,y[X) ~ A^/(2T(t^+l)) such that the integrand diverges in this limit. Therefore, we 
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res6 



res7 



res8 



include an additional factor k^^^^f^ in the logarithm in order to avoid contributions 
coming from the origin. Because there is no additional pole enclosed, this does not 
change the result. Furthermore we should note that the integrand has branch cuts 
starting at A = ±«(p^ + m?). Leaving out the n, £ summations for the moment and 
considering the A-integration alone, we then obtain {ly ^ £ + ^) 



CPMi^) 



^^ipl + X^ + mT'^HX-^J.W] 



(3.8) 



7 

sinTTS 



n 



d\ 



dk{e-pl-mT^^\n[k-^'+ih,^^{k) 



VpI+^' 



where J^{ik) = e*^''J^(-zA:) and I^{k) = e-"'''/^J^{ik) has been used fIE\ . 

The next step is to add and subtract the asymptotic behavior of the integrand in 
p.8p . The relevant uniform asymptotics, after substituting k = vz in the integral, 
is the Debye expansion of the Bessel functions [1] . We have 



(3.9) 



Iy{vz) 



1 



'2wu {I + z2)4 



E 

fe=l 



Uk{t) 



with t = l/vTT~2^ and rj = Vl + 2^ + ln[z/(l + Vl + z'^)]- The first few coefficients 
are listed in [1] , higher coefficients are immediately obtained by using the recursion 

m 

t 

(3.10) Uk+l{t) = ^t^{l-t^)u[it) + ^ fdT (l-5T2)ufc(T), 



res9 



reslO 



starting with wo(i) = 1- As is clear, all the Ufc(t) are polynomials in t. The same 
holds for the coefficients D„(i) defined by 



(3.11) 



In 



1 



E°° Uk(t) 



E 



Dn{t) 



The polynomials Ufc(<) as well as Dn(t) are easily found with the help of a simple 
computer program. As we will see below, we need the first three terms in the 
expansion (J3.11I) . Explicitly we have 



(3.12) 



Diit) - 
D2it) - 



1 5 . 

8 24 ' 

-t^ - -t^ 
16 8 



25 



:f 



531 



16 



221 



-t' 



1105 



-t" 



384 640 128 1152 

Adding and subtracting these terms in p.8p allows us to rewrite the zeta function 



as 



Cp(s) =CP,/(s)+CP,a.(s), 
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reslOa 



reslOb 



resll 



reslla 



resllb 



where 



oo oo 






Pn-'m- 



2\-s 



n— — oo (t—Q 



dz 



(3.13) y.— {\n[z~''I^{vz)]^\n 



Di(t) D2{t) Dsit) 



Cp.asis) 



(3.14) 



2^(l + z2)l/4_ 
oo 

Y^ Y,{2£+1) f dz{z^v^-pl-m^)- 



OO OO 



n— — oo ^—0 



dz 



/2^ (1 + ^2)1/4 



Di(f) , i?2(i) , D:i{t) 



The number of terms subtracted in (j3.13p is chosen such that C,pj{s) is analytic 
about s = 0. The contributions from the asymptotics collected in (j3.14p are simple 
enough for an analytical continuation to be found. Although it would be possible 
to proceed just with the contribution from inside the ball, in order to make the 
calculation as transparent and unambiguous as possible (as far as the interpretation 
of results goes) let us add the contribution from outside the ball. 

The exterior of the ball, once the free Minkowski space contribution is subtracted, 
yields the starting point p.Sp with the replacement k"'^!^, — > k'^ K^, [5]. In this case 
the relevant uniform asymptotics is [T] 



(3.15) 



Ky{vz) 



-vrj 



2i^{l + z^)i 



i + E(-i) 



',Ukit) 



fc=l 



where the notation is as in p.9p . This produces the analogous splitting of the 
zeta function for the exterior space. Due to the characteristic sign changes in the 
asymptotics of I^ and K^, , adding up the interior and exterior contributions several 
cancelations take place. As a result, the zeta function for the total space has the 
form 



Ctot(s) = Ctotjis) + Ctot,asis) 



with 

Ctotj{s) 

(3.16) 

Ctot,as{s) 

(3.17) 



OO OO 



E E(2^+i) / 



1/22 2 2n- 

dz [z V — p^ ~ m ) 



n= — oo 1=0 



^ |ln[/,(i/z)A',(i/z)] + ln(2i^) + i ln(l + ^') - ^ 



D2it) , 



OO OO 



E E(2^+i) 



n— — OO i—0 



d_ 

dz 



ln(2i/) ln(l + z 



J / 2 2 2 2\' 

dz [z V — p^ ~ m ) 



r^2(i) . 
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By construction, Ctot,/(s) is analytic about s = and one immediately finds 



resl2 



Ctotjio) 

(3.18) 



OQ OO 



E E(2^+i) 



n— — OO i—0 



In (/,(j/0)X,(i/z)) + ln(2i/) + i In (l + z^) - ^D^it) 



with t — 1/Vl + z^ as defined earlier. Although one could use p.lSp for numer- 
ical evaluation, further simplifications are possible. Following [14| we rewrite this 
expression according to 



resl2a 



(3.19) 



1 + ^2 ^ 1 + Pzi. 



1 + % 



The advantage of the right hand side is that it can be expanded further for v^ -^ co 
or p^ — )■ (X) or both. This will allow us to subtract exactly the behavior that makes 
the double series convergent; the oversubtraction immanent in p.lSp can then be 
avoided. It is expected that expanding the rightmost factor further for i/^ +p^ ^ 1 
leads to considerable cancelations when combined with Ctot,as(0) El- 

We split the asymptotic terms in (|3.18p into those strictly needed to make the 
sums convergent and those that ultimately will not contribute. For example, we 
expand according to 



ln(l + z2)| 



In 1 



In 1 



vl 



In 1 



In 1 



vi + m^ 



1/2 + p2 



V + Vn 



In 1 



v'+vi 



Vn 



The first two terms have to be subtracted in ()3.18p in order to make the summations 
convergent. The terms in square brackets are of the order OiXjiiP' +p^)^) and 
even after performing the summations in p.lSp a finite result follows. Thus the 
first two terms represent a minimal set of terms to be subtracted in p.l8|) in 
order to make the sums finite. This minimal set of necessary terms will be called 
In/^^^™ (z-y/p2 _|_ 777,2 )_ ^\\fi last two terms can be summed separately yielding 
a finite answer; those terms are summarized under In/"''^'"'^ (*\/Pn + 'm?'). One 
can proceed along the same lines for all other terms. With the definition 



In Z^"'*^™ (i Vrc + m2 ) 



In(2i.)-lln(l + z2) + ^i52(i) 
2 v^ 



resl3 (3.20) 



ln/r™-(^)(zv/^^TT;?) + ln/r™-(^)(zv^^T^) 
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11 



the splitting is 



in/r-(^'(*v^^T^) = -i-(2-)-^i^^(i + §)-^^;2^ 



resl4 



(3.21) 



1 1 



1 



16 



ln/r^"'(2)(zV5^T7^) 



Pl\ 8A^p|A2 16 
1 m^ 



P„ 



1 / m^ 



2 1 1 



v^ 16 



1 



^ '^Pn 



resl5 



(3.22) 



1 \2 



v^+f 



16 ^ + C^ 



u^+Pi 



We have used the given notation for the asymptotics to make a comparison with 
[l4] as easy as possible. With these asymptotic quantities we rewrite Q^^ f (0) as 



oo oo 



resl6 



n— — oo i—0 

X [ln(/,(Vp2+m2)i^,(Vp^ + m2)) - ln/r^'"'«(z Vp,U m^)' 

OO OO 

(3.23) + Yl E(2^+l)l^/r"'^'^(*V^lT^)- 



n— — oo ^^0 



Let us next analyze C'tot asi^)- To analyze Ctot,as{s), equation p.l7p . further we use 
the integrals 



du \ u^ I sm TTs 



oo 

/ .„(.=-p;.-.»=)-^(i4 



2\ -N/2 



vpI+^ 



7rr(s + 



N^ 



-2s 



sin(7rs)r (f ) r(s) / ^ p^^y 
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resl7 



resl7a 



reslS 



resl9 



which are the relevant ones after substituting zv = u. This shows 



C30 OO 



c*o,a.w = - E E-'^'Mi 



p„ + m^ 



2 I 2\ — s— 1 



(3.24) 



n— — OO i—{) 

OO OO 

n=— OO £=0 

-^^ E E^-^-^-^i 

n— — OO ^^0 
5S(S + 1)(S+^ y. y. 



OO OO 



8 



Vi + JTi^ 



2 I 2\ — S--3 



ri=— OO £=0 



To each of these terms we apply the rewriting p.l9p . Intermediate expressions are 
relatively lengthy and we explain details only for the first term. We proceed as for 
the splitting in p.2ip and (|3.22p and write 



OO OO / 9 , 9 \ — S 

2 \ -s 



-2s / -. . Vn 
9 



-E E-^ 

n— — OO £—0 

OO OO 

= - E E-^ 

n— — OO £—0 

OO OO 

= " E E (;y2+p2)s 

OO OO 

= - E E 



1/2 +p2 



n— — OO t—^ 

OO OO 



(^2+p2) 



,2^s 



(3.25) -i E E 



2v 



2 ^ f-(j^'+P^)" 

n=-oo £=0 ^ -^"' 



STTi 



m2 y 


m2 
1 1 " 


^^+pj 


' 2 , 2 


m2 \-^ 


11" "^' 


i.2+p2J 


'"-^+P?. 


2 OO OO 


2zy 



1-s 



1/2 +p2 



n— — OO ^^0 



(!/2+p2).+ l 



Note, that the first line is seen to be analytic about s = 0. We have subtracted 
the minimal number of terms to make the sums convergent. The remaining terms 
represent Epstein type zcta functions, 



OO OO 



(3.26) 



i;W(.,a)= E E2^ 



n— — OO £^0 



(i/2 + a^v?')^ ' 



the analytical continuation of which is well understood. Performing a Poisson 
resummation in the n-summation [31 1151 130] yields 



i?W(s,, 



20Fr(s-i) 



(3.27) 



a r(s 



Cff(2s-fc-2,l/2) 



^ ' £=0 



mL^'^-'Y.^-'^H-si:^^ 



The first line has poles at s = 1/2 — j, j E Nq, and for s = (fc + 3)/2, the second 
line is analytic for s £ C 
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res20 



In terms of these Epstein functions, in equation p.25|) we have shown 



oo oo 



(3.28) 



OO CX3 



.eW 



!,■ 



27r 



E' 



,2 As 






Pn 



-1 



P« 



Noting from equation p.27p that E^^'>{s,a) and E^'^^s + l,a) are analytic about 
s = 0, we get 



oo oo 



res21 



d^ ^ ^^ 

n— — oo i—0 



(3.29) 



Y-j,i-2s 



1 



P« 



OO oo 



^ 4-'°"(°.f 



In 1 



2 V /3 



Pn 



ly^ + pI 



res22 



The last term on the right cancels the first line from In / 



asy7n,(2) 



[iy/pi + Tri?) in 



equation p.22p . the remaining terms are easily found from p.27p . 

One can proceed in exactly this way for the other terms in Ctot.as(s); there 
are always terms that cancel with terms from In/^'*^™ (i-^/p^+rT?) in equation 
p.22p and terms expressible using the Epstein type zeta functions given in equation 
(|3.26p . Adding up all contributions, the second line in equation p.23p completely 
cancels and we obtain the following closed form for the finite-T zeta function, 



CL(o) 



oo oo 






Y,2v [ln/,(v/p2+m2)i^,(v^ 



8 i^2 + p2 



Pn 



^v'^ 



^..o.,(„,-),l„,,..,(,|) 



(3.30) 



4 



' /?; 8 V /?; 



From equation p.27p it is clear that the Epstein type zeta functions contain zero 
temperature contributions to the Casimir energy (first line in p.27p ) and exponen- 
tially damped contributions for small temperature described by the Bessel functions 
(second line in p.27p ). As it turns out, the zero temperature contributions from the 
Epstein type zeta functions in p.30p all vanish. The remaining zero temperature 
contributions in (|3.30p are found replacing the Riemann sum over n by an integral. 



oo °° 

E /w =^ 1^ / rfp/(p)- 
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As /? — >■ this shows 



«=0 



res23 | (3.31) + ln{v^ + p'^) + 



from which the Casimir energy 
simpler than previous results given 
be performed. 



]2z^ In/, 




iVp' 


+ m' 


)if.(v^ 


2 + 7^2) 




m2 


1 


1 


1 


61/2 

1 


5i/4 Y 


^2 _j_p2 


81/2 


+ p2 


I (^2 


+ p2)2 


l/2+p2)3JJ 



is trivially obtained. The result is much 
1 [5] and a numerical evaluation could easily 



res24 



res24a 



res25 



4. Functional determinant on a two dimensional torus 

As our next example let us consider a two dimensional torus S^ x 5*^. For 
convenience we choose the perimeter of the circles to be 1. The relevant eigenvalue 
problem to be considered then is 

92 ^2 

and we choose periodic boundary conditions 

(l>\{x, y) = 0(a; + 1, y), (t)\{x, y) = (j)x{x, y + 1), 

d(j)x{x,y) d(j){x + l,y) d(j)x{x,y) d(j){x,y + l) 



P4'\ix,y) 



(l)\{x,y) = A^(/)a(x,2/), 



dx 



dx 



dy 



dy 



Eigenfunctions and eigenvalues clearly are 



4'm,n{x,y) = e 



-27Timx „ — 27rmv 



A^ = (27r)2(m2 + n2), n,me 



The related zeta function then reads 

(4.1) Cp(s) - (2^)-^-' 



y^ (to2 + ^2)" 

(m,n)eZ2/{(0,0)} 



note that the zero mode m = n = has to be omitted in the summation to make 
Cpis) well defined. The zeta function in equation (j4.ip is an Epstein zcta function 
and Cp(0) can be evaluated using the Kronecker limit formula [18l[T9]. Here, we 
apply the contour approach previously outlined which simplifies the calculation. 

Instead of using the fact that the eigenvalues can be given in closed form, we 
proceed differently. We say that 

A^ = (27r)2(n2 + fc2)^ n e Z, 

where k is determined as a solution to the equation 

(4.2) e"''^ - 6""'= = 0. 



Of course, solutions are given by k € Z and the correct eigenvalues follow. Using 
equation ()4.2p determining the eigenvalues in the way we have used equations (|2.3p 
and p.4p . the zeta function can be represented as the contour integral 



Cp(s) 



4E 



dk 



{2n)-^'{n' + k^y 



, d , / e 



■nik 



-Txik 



2iTik 



(4.3) 



+4(27r)-2^Cfl(2s). 
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The second line represents the part where one of the two indices tti or n is zero in 
equation (|4.1[) . The first fine represents the remaining contributions. The factor of 
4 is a rcsuh of summing over positive n only and because the contour 7 is supposed 
to enclose positive integers only. The reason that we have used 

pirik ^ — Trik 



2Trik 
instead of the above equation (|4.2p is that 

p-Kik p — Txik 

lim = 1, 

fc-yo 2TTik 

which will allow us to shift the contour in a way as to include the origin; see the 
discussion below equation (j3.7p . Let us evaluate the contour integral 

r rlh rl / p-rrik _ -Trik 



27ri dk \ 2Tiik 

7 



Substituting k = ^/z and deforming the contour to the negative real axis along the 
lines described previously, an intermediate result is 



00 



Tii26l (4.4) C„(s) = (27r)-2^^^^^ /dz(z-n2)-^4ln 

TT J dz 



d , /e'ryi_g-7ryi 



2Try/z 



From the behavior of the integrand as z -> 00 and z — >■ n^ this representation is 
seen to be valid for 1/2 < Sfts < 1. In order to construct the analytical continuation 
to a neighborhood of ,s = wc note that 






2'K'Jz 21x^/2 

We therefore write 

_2sSin7rs 



C„(s) = (2T,y \ dz{z-n')-'—\M 



d , I e^^ 



dz \ 2'K'Jz 



_^^^y2s^2!lll I dziz^nT'^^Jn[l-e-''^^ 



The first line is evaluated using 



dz^ — ^ = ^^r(i-s)r -- + S , 



, (z-n^)-' 
dz 



With the identity [28] 

sinTTS 



SmTTS 



r(i-.) 



r{s) 
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this produces 

Cn(s) 



,ni-^^r(-i + .) 1 



(27r)-2'*n- 



l(27r)-'"-^ 

oo 
TT J dz \ 



This is the form that allows the sum over n to be (partly) performed and it shows 



Cp(s) = 4(2^)-2^5^^^y dz{z -n^)-'^ In (l~e-^^^ 
n ^-^ I dz \ 

n—l 9 



rees27 



(4.5) 



0Fr(s) 



~2(27r)-2-Cfl(2.s) + (2^)- 



&(2s-l). 



This form allows for the evaluation of Cp(0)- From known elementary properties of 
the F-function and the zcta function of Riemann 1281 we obtain 



d 
ds 



[{2it)-^%r{2s)\ = -2 ln(27r)Cp(0) + 2(^(0)^0, 



d_ 

ds 



(2^)_2.+i£Li±il^^(2.-l) 
V^F(s) 



TT 

3' 



The first line in (|4.5p is also easily evaluated because 

oo 
oo „ , 

, ,„ , 9„sm7rs v^ / , / 9n « " , / 

4(2^)"'^-— E dz{z^nY'-\n[l^e 
=0 ^ „=i -4 "^^ 

CX3 „ , OO 



res28 



n—l o ?i— 1 

This can be reexpressed using the Dedekind eta function 

oo 

77(t) :== e"^/i2 [] (1 _ e^""^) 

n=l 

for T e C, 3fiT > 0. The relation relevant for us follows by setting t = i, 

oo 

lnh(0|* = -|+4^1n(l-e-2™). 

^ r.=l 

Adding up all contributions for Cp(0)i the final answer reads 

oo 

(4.6) Cp(0) - f - 4^] In (1 - e-2"") = - In |r;(z)|4 

■^ n=l 

in agreement with known answers; see, e.g., [iSllSO] . 
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5. Conclusions 

In this article wc have shown that contour integrals arc very useful and effective 
tools for the evaluation of determinants of differential operators. Although the 
results look very simple only in one dimension, see equation (|2.10[) . for particular 
configurations also in higher dimensions closed answers can be found suitable for 
numerical evaluation, see equations p.3ip and (|4.6|) . Here we have provided answers 
only for the torus and a spherically symmetric situation. But the same ideas should 
apply when separability of the partial differential equations in other coordinate 
systems is possible. Results in this direction will be presented elsewhere. 
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